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Riassunto. In questo seminario presentiamo alcuni risultati ottenuti in collaborazione con 
M.C. Tesi sulla costruzione di uno schema agli elementi finiti per il problema di Dirichlet 
in un rettangolo del piano asssociato a un operatore ellittico degenere ‘tipo Grushin”. Il 
metodo consiste essenzialmente nell’utilizzare una triangolazione che rispetta la geometria 
associata all'operatore tramite la metrica di Carnot-Carathéodory. 


Abstract. In this seminar we sketch some recent results obtained jointly with M.C. Tesi 
concerning a finite element method for the Dirichlet problem associated with a degenerate 
elliptic operator of Grushin type. Basically, the method consists of using a triangularization 
fitting the Carnot-Carathéodory metric associated with the operator. 
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I risultati esposti in questo seminario sono stati ottenuti in collaborazione con M. C. Tesi 
([FT]). Sia Q=]-1,1[x]-1,1[uncuboinR?=R,xR,, e sia il suo bordo. Consideriamo 


l’operatore del secondo ordine in forma di divergenza in Q definito da 


3°, 
(1.1) L=-), d(a;(2)9;), 


1,j=1 
dove i coefficienti a;; = a;; sono funzioni reali misurabili e, per qualche v € (0, 1), 


2 
(12) VE + (299) < DO assai < TE + XP). 
tu=i 
per ogni 6 = (£,n) e 2 = (2,y) € R?. Qui A è una funzione Lipschitziana nonnegativa 
in R. Operatori di questa forma sono noti in letteratura come operatori tipo Grushin e 
le proprietà di regolarità delle soluzioni deboli di Lu = f sono state ampiamente studiate 
negli ultimi anni: si vedano, ad esempio, [FL], [X], [FS],[F .[FGuW1,2]. 

Un produttivo approccio allo studio di questi operatori ha consistito (si veda [FL]) 
nell’associare all’operatore £ una opportuna metrica (non Riemanniana) d che è fonda- 
mentalmente data dal minimo tempo necessario per passare da un punto a un altro lungo 
curve continue che sono a tratti curve integrali dei campi vettoriali +0; e +0, (si veda la 
Definizione 2.1 per una definizione precisa). Se per esmpio siamo interessati alla regolarità 
Hélderiana delle soluzioni deboli (teorema di De Giorgi-Nash-Moser) o alla disuguaglianza 
di Harnack per le soluzioni deboli positive, possiamo ripetere le tecniche classiche svilup- 

‘ pate per equazioni ellittiche ([DG], [Mo]) sostituendo le usuali sfere Euclidee con le cosidette 
sfere metriche, cioè le sfere della metrica d. À 

Lo scopo di questo lavoro è quello di mostrare che un simile approccio geometrico 
ci permette di sviluppare una approssimazione con elementi finiti per questa classe di 
operatori. Infatti costruiremo una triangolazione di una regione piana per mezzo di una 
famiglia di triangoli non isotropi che rispecchia la geometria associata alla metrica d, nel 
senso che ogni triangolo della nostra triangolazione contiene ed è contenuto in due sfere 
metriche di raggi confrontabili. La forma di questi triangoli non sarà banale a descriversi, 
poichè le sfere metriche non sono invarianti rispetto alle traslazioni euclidee, così che non 
possiamo semplicemente ripetere una stessa sfera traslandola e, analogamente, non ci sono 
omotetie che permettano di riscalare le nostre stime. 

Recentemente, nello stesso spirito un metodo alle differenze finite per equazioni tipo 
Kolmogorov è stato sviluppato in [MP]. 

Sottolineamo che il nostro approccio non è precisamente un metodo adattivo, perchè, per 
così dire, la geometria è fissata a priori ed è data dal nostro operatore modello 0? + \2(x)d? 
che gioca il ruolo di un operatore di Laplace-Beltrami per la nostra geometria. Anzi, un 
metodo adattivo può essere sovrapposto alla nostra geometria, in modo da tener conto 
dell’oscillazione dei coefficienti. 

Notiamo esplicitamente che, a causa della mancanza di ellitticità dove A si annulla, dob- 
biamo cercare soluzioni deboli in spazi funzionali più grandi dell’usuale spazio di Sobolev 


É 1(9) e che sono dati dal completamento di C°(0) rispetto alla norma 
ulizzca + N u]lz2) + A92u]|z20); 
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così che in generale le nostre soluzioni deboli non appartengono a H!(9Q). 


Questo approccio è stato usato per una classe molto più grande di operatori ellittici 
degeneri, il cui prototipo è dato dagli operatori somma di quadrati di Hòrmander in R" 
della forma Vici Xi, dove X1,..., Xn sono campi vettoriali infinitamente derivabili tali 
che il rango dell’algebra di Lie generata da essi è uguale a n in ogni punto. Ad esempio, 
se scegliamo Me) = |e|F, per qualche naturale &, il nostro operatore modello 0? + |e|?F02 
è un'operatore di Hérmander. Poichè noi tratteremo con funzioni ) non infinitamente 
derivabili, dovremo imporre ulteriori condizioni su À per sostituire questa ipotesi del rango 
(si veda l’Ipotesi (H) più sotto). 


Se noi cerchiamo di seguire lo schema della prova di Moser della regolarità puntuale delle 
soluzioni deboli di Lu = f, fin dal principio ci sono due punti che appaiono giocare un ruolo 
chiave: il fatto che la metrica d è doubling (cioè il volume della sfera metrica di raggio 2r è 
controllato da una costante per il volume di una sfera metrica di raggio r con lo stesso cen- 
tro), e una opportuna disuguaglianza di Sobolev-Poincaré sulle sfere metriche, dove a destra 
dobbiamo sostituire l’usuale gradiente Vu con il ‘gradiente degenere’ Viu = (01,403) as- 
sociato all'operatore. Queste disuguaglianze contengono profonde informazioni riguardanti 
la geometria associata alla metrica d, poichè mostrano che la dimensione geometrica dello 
spazio metrico definito da d è molto più grande di 2 (o di n in generale) e, per così dire, è 
tanto più grande quanto più \ è degenere. Questo fenomeno è stato studiato nel contesto 
più generale dei campi vettoriali tipo Hòrmander e appare chiaramente in una famiglia di 


disuguaglianze isoperimetriche (si vedano [FGaW], [FLW], [CDG1,2], [GN], [Gr]). 


Purtroppo questo fenomeno dimensionale influisce negativamente sulle nostre stime 
dell’errore. Infatti, per prima-cosa, non abbiamo un teorema di immersione di Sobolev 
per controllare i valori puntuali di una soluzione debole per mezzo di norme di Sobolev 
di ordine s > 1 = n/2, per n = 2, come nel caso ellittico, proprio perchè la dimensione 
di (9, d) è in generale molto più grande di 2. È tuttavia possibile superare questa diffi- 
coltà, ma lo stesso fenomeno dimensionale riappare nell’approssimazione numerica, perchè, 
in corrispondenza a una griglia di n? punti, troviamo una stima dell’errore dell’ordine di 
n1/(1+1), dove y > 0 e y+2 è la dimensione geometrica di (9, d). In altre parole, ci occorre 
un grande numero di triangoli per ottenere piccoli errori, molti di più che nel caso ellittico, 
e tanti più quanto più A è ‘piatto’ nei punti dove si annulla, così che la nostra approssi- 
mazione converge, ma l’ordine di convergenza dipende dall’ordine di degenerazione di À. 
E quindi necessario tener conto di questo fenomeno quando si confrontano i nostri risul- 
tati numerici con quelli che possiamo ottenere semplicemente facendo girare ‘brutalmente’ 
‘procedure numeriche fuori da ogni schema teorico ottenendo buoni risultati fuori dagli zeri 
di \ (poichè l'operatore £ è localmente ellittico in queste regioni). Discuteremo anche 
l’implementazione algoritmica del metodo, e mostreremo, analizzando opportuni esempi, 
che — come è giusto attenderci - la stima dell’errore nella norma dell’energia può essere 
migliore di quella ottenuta con una triangolazione uniforme, o anche con una triangolazione 
uniforme adattiva (ma sottolineiamo ancora che l’uso di una triangolazione uniforme non 
ha giustificazione teorica, perchè ci sono soluzioni che non appartengono agli usuali spazi 


(c) 
di Sobolev H!(Q)). Per finire, forniremo esempi numerici che mostrano che la nostra stima 
dell’errore è ottimale. 
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2. In questa lavoro denoteremo un generico punto di R? con z = (x,y). Nel seguito 
assumeremo che la funzione A soddisfi la seguente ipotesi: 


Ipotesi (H) Esiste una costante positiva c1 tale che, per ogni intervallo compatto I CR 


1 
KA <- & 
0<ci maxÀ & mf S maxA 


dove |I| è la lunghezza di I. 


Questa condizione è chiamata condizione RH > in [F] e [FGuW1]; per commenti sul suo 
significato geometrico, si veda [CF]. 


Definizione 2.1. Diremo che una curva assolutamente continua Y: [0,7] + R? è una 
curva sub-unitaria se per ogni ( = (€&,n) € R?, 


< It), >?< JE? + A°(1(8)Im]? 


per q.0.t € [0, T] (per semplificare le nostre notazioni abbiamo considerato qui À come una 
funzione di 2 € R?). Se 21,2, € R?, poniamo 


d(z1,22) = inf {T > 0; esiste una curva sub-unitaria 7 : [0,T] + R? 
tale che y(0) = z1,9(T)= 22}. 


Per l’ipotesi (H), d(z1,22) < 00 per ogni 21, 22 € R? e quindi è una metrica. 
Se ora z = (x,y) e R? er > 0, poniamo 


t+ 


(2.1) F(z,r)= F(2,r) -J 9ds 


T 
Diremo che una costante c > 0 è una costante geometrica se dipende dalla costante ci 
dell’Ipotesi (H) e da sup A. Per alleggerire le nostre notazioni, talora indicheremo con la 
stessa lettera c differenti costanti geometriche. 
Abbiamo: 


Proposizione 2.2 ([FGuW1], Proposizione 2.5). Supponiamo che valga (H). Allora 
per ogni punto zo = (r0,yo) € R? esistono un intorno U di zo e una costante geometrica 
Y2 0 tali che 
(i) F(2,0t) > c@!+1F(2,t),0<0<1; 
(ii) F(2,0t) < c@+1F(z,t),0>]1; 
(iii) F1(2,0t) < c01/0+VF-1(2,t),0<0<1; 
(iv) F1(2,04) > c@/0+) F-1(2,t),1<0; 
(v) F(2,0t) < c0F(z,t),0<0<1; 
(vi) F(2,0t) > c0F(z,t),0>]1; 
(vii) F_1(2,0t) > c@F-!(z,t),0<@< 1;in particolare, se r appartiene al compatto 
|1,1), allora F(z,r)<cr perr < f; 
(vii) F71(2,0t) < c0F!(2,t),1<0; 
(ix) se d(z1,z2) < ct, alloraF(z,,t) > F(22,1), 
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per 0<t<t e z EU, dovec è una costante geometrica. 
In particolare, da (i) segue che 


(2.2) îI Mx + sé) ds > ct1+? 


per ogni <= (2,Y) È U, É E So = {£ eR ; E — fol S 6} Cc [-1,1]\ {0} ete (0, to). 
Osserviamo che dalla Proposizione 2.2 (iii), si ottiene la seguente proprietà di doubling 
(2.3) F(z,r) < F(2,2r) < cF(z,r) 
per ogni 2 E YU, r < To, ro e c essendo costanti geometriche. 
Teorema 2.3. Supponiamo (H) soddisfatta. Allora: 
(i) d(21,22) < 00 per ogni z1, 22 € R° e quindi d è una metrica. 
(ii) se denotiamo con B(z,r) la d-sfera di centro 2 e raggio r (cioè B(z,r) = {2' € 
R?; d(z,2')< r}), allora esistono due costanti geometriche t, > 0 e b > 1 tali che, 
per ogni zo € R? e t € (0,t1), abbiamo: 


(2.4) Q(z0,t/0) C B(z0,t) C Q(zo, dit), 
dove, per ogni r > 0, 
Q(z0,7) = {2= (2,9) e R°;le — col <r e ly— yo < F(z0,7)}. 
(iii) esistono due costanti geometriche A > 0 e ro > 0 tali che 
(2.5) |B(,2r)] < A|B(G,)] 


per ogni 2 € U er € (0,ro), cioè lo spazio metrico (R?, d) è uno spazio di tipo 
omogeneo rispetto alla misura di Lebesgue. 


(iv) Se@>0, allora 
(2.6) c1(9)|B(2,r)] < |B(2,9r)] < c2(9)|B(z,r)] 


per ogni z € U e r < ro e per opportune costanti c1(9) e c2(9) che sono costanti 
geometriche,eccetto per la dipendenza da 0. 


Nel seguito, denoteremo con Vy = (4, Ad») il gradiente degenere associato all’operatore 
L e porremo 


(2.7) |Vau]? = [Gru]? + A°(2)|Q2u]?. 
Inoltre, denoteremo con H1(Q) lo spazio di Sobolev degenere associato con Vi, cioè 
l’insieme delle u € I?(9) tali che 


du € L°(Q) X0zu € L?(Q), 
fornito della norma naturale 
(2.8) lullin = luz: + NA]l7(0) + [AG2u]{7:(0) 
Notiamo che in questi spazi vale un teorema tipo Meyers-Serrin, cioè 
C®”(9) N Hi(9)è denso in H}(9) 


(si vedano [Fr], [FSSC] and [GN]). Quindi sarà naturale denotare con H1(0) la chiusura 
di CE(Q) in Hi(2). 
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Teorema 2.4. Supponiamo valga (H); allora esiste una costante geometrica c > 0 tale 


che 
fi lu|?dz < cf |Vau|?dz 
Q (e) 


Possiamo ora formulare il nostro risultato principale riguardante il problema di Dirichlet 
per Lin Q. 


Teorema 2.5. Siano fo, f = (fi, f2) tali che fo,|f| € L?(9). Allora esiste una ed una 
sola u € H}(9) che risolve il problema di Dirichlet 


per ogni u € H1(Q). 


Lu= —di in Q 
(2.10) (2) { sessi 
u=0 su T, 
dove divyf = di fi + A@; fa, nel senso che 


A(u,9)= É {991 9+ X°dyudyp}dz = È {f1016 + fa0x.9 + fog}dz = Ly(0) 


per ogni p E CE°(Q). 
Ragionando come in [F] e [FS], Teoremi 5.11 e 6.4 rispettivamente, possiamo provare il 
risultato seguente. i 


Teorema 2.6. Se fo, |f| € L? per p > po = Po(A), allora la soluzione u di (P) è Holderiana 
inQ. 


3. Possiamo ora descrivere la nostra triangolazione. 


Teorema 3.1. Per ogni n > 0) esiste una decomposizione finita Tr del dominio © 


dove 


(i) ogni K è un triangolo compatto con interno non vuoto Int(K); 
(i) Int(K1)N Int(K2)=0 perK,,K2€T, distinti; 
(ili) se F = K,N K. # 0, dove K e K, sono elementi distinti di T,, allora F_ deve 
essere un lato per entrambi K, e Kz o un vertice per entrambi K, e Ka; 
(iv) per ogni K € T, esistono zx, îK € R?, Fx, FI > 0 tali che 


B(zx,Tx) CK E B(Ex,Fk), 0<e< È <c, 


c e C' essendo costanti geometriche; 
supy Fr < const n-1/(1+%), dove Y è la costante della Proposizione 2.2. 


(v 


n 


90 


Cenno della dimostrazione. Limitiamoci a costruire i vertici della triangolazione in Q 
dove x > 0,y > 0. I rimanenti vertici si otterranno per riflessione. 
. 1 È 
Scegliamo a > 0 tale che a fn A(s)ds = 1. Possiamo sempre assumere che la costante 


ci in (H) sia tale che aci < 1, e siano 60, 61, -.-, n scelti in modo che 
Siti 1 
(3.1) 60 = 0 af A(s)ds= = for j=0,1,...,n-1. 
5; a 


A tal fine basta porre A(t) = A A(s)ds ( che è strettamente crescente per (H)), e 


(3.2) 6;= A (1), fel 


an 


Consideriamo allora la triangolazione di Q* data dalla famiglia di nodi 


Notiamo che i nodi 


stanno su T. 


Possiamo ora procedere in modo standard per definire una spazio di funzioni di di- 
mensione finita V,, nel modo seguente: siano P,,...,Py,N = N(n)i nodi di 7 che 
appartengono a Int(2). Consideriamo l'insieme ®, di tutte le funzioni continue lineari a 
tratti 9;,j=1,..., N tali che 


(3.5) g;=0onT and pg;(Pi)=d; i=1,...,N, 
e denotiamo con V, lo spazio lineare generato da Dn. 
Lemma 3.2. VW C Hi(Q) pern> 1 
Dimostrazione. Basta notare che V, C H1(9) C HX(0). 
Una funzione vn € V ha ora la rappresentazione 
N 
(3.6) vn(2) = Dl vn(Pi)pi(2) 
i=l 


e il nostro problema di Dirichlet può ora essere approssimato dal problema seguente: 


91 


trovare un E V,, tale che 
(3.7) (Pi) A(Un;0n) = Ly(0n) Von E Wa 


Come nella teoria ellittica, questa problema può essere risolto risolvendo un sistema lineare 
N x N la cui matrice di stiffness ha elementi dati da (AP: 9;)); ja SP 


Ancora come nella teoria classica, possiamo definire un operatore di interpolazione II, : 
CA) —+ Va ponendo 


N 


(3.8) i In(v) =) uPi)g: 


i=1 


4. Enunciamo la nostra stima dell'errore sotto ipotesi un po’ più restrittive su \. 
Teorema 4.1. Supponiamo che \ sia una funzione C®1 che soddisfa (H) e tale che X? € 
Cl, i suoi zeri siano isolati e appartengano Int(Q), e che (A?)" > 0 in un intorno di questi 
zeri. Se f € L°(Q), allora esiste una ed una sola u € H1(9) tale che 

Lou=f inQ. 


o 
Inoltre, le approssimazioni di Galerkin un € V, definite da (3.7) convergono au € H 19) 
per n + co e vale la seguente stima dell'errore: 


lu unllaza) £ Bo-VO+A FIL: 
dove B è una costante geometrica. 


Dimostrazione. Possiamo sempre supporre che M(0) = 0 e Mx) > 0ifr #0. Inoltre, 
possiamo assumere che M(z') > M(r") per 0 < 2" < c'e Me) < Mx") pere" <a <0. 
Come in [QV], Teorema 5.2.1, esistono B,, B. > 0 dipendenti solo da v e dalla costante c 
di (2.9) tali che 


(4.1) urla) < Billfilcz() 
2 
(4.2) un — ul) < B2 siRE, Ilu— vallata) 
D’altra parte 
(4.3) n" lu valli) £ lu — In()]7:(0) 


poichè IIn(u) è ben definito in forza della continuità di u e appartiene a V, per il Lemma 
3.2. Combinando (4.2) e (4.3) segue che un + u in H}(9) una volta provato che lu — 
In(u)] #7!) —+ 0 pern4 co. 

Ora 


lu En(u)Il aio) S c{Ilu—Tn(u)]lz2(2) + |Va(u — Tn(u))]z2(0)} 


(4.4) =c ) {llu- Un(u)llz:(x) + [Va(u— In(u))]z2)} 
KET, 


Abbiamo: 
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Lemma 4.2. Esiste una costante geometrica B3 > 0 tale che, se poniamo X1 = di, 
X,= A@, allora, se K € Tn 
lu Tn (4)]lz21) + |Va(u — In(u)llz2 0) £ 
(4.5) B3 {ullz20) + |XZullezo + 1X2X1ellz:4 pereeta, 
Notiamo che Int(K) non interseca nè I nè la linea di degenerazione r = 0, così che 
u € C?:"(K), per le stime classiche di Schauder. 


Proof. Prima stimiamo 


3 
46) f [Wa(u- tn(u)ldzdy = f,Iva@- Du0)eld = de, 


i=l 


dove Q1, 02,23 sono i vertici di K. 
Supponiamo ad esempio che che K abbia vertici (8; E), (641) E), (6;, EL), e poniamo 
r=68+(5;41—- G)r 
(4.7) { _ : ( i i) 
vertal- 


Questo cambio di variabili muta K nel triangolo fisso di vertici (0,0), (1,0), (0,1) e il suo 
determinante jacobiano vale 2|I], così che, se denotiamo con è una generica funzione v 
scritta nelle nuove variabili x',y'", otteniamo 


3 
(4.8) f.to tinto) Paedy = 21k1 f,18- Donata 
K K i=1 


Notiamo ora che i 
O,:0 = (6j+1— 6;)0zv Oy 6 = > Bon 


Dunque 
3 T_T 
f, jax(u— YO u(Q:)pi)dedy = 2] SI i) 
L i=l 
2|K] 3% dt 
=T—_T- 0z'(...)|fdr' dy 
(841 — 5;)? A (o 
° 1 i 7 292!4,) 
4.9 = tr 0: (...)|fdr dy . 
va n(6;+1 — 5;) Ai (H 
Analogamente 


3 


i DIL — YO Qaeda 


= 
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= (fia = 85)  A°(8;+ (6541 — Gay da' ay 

K 
< (6541 — 65): (max A f lan laz'y 

(5; ,5; +1] K 

1 î 6; 41 2 SE 

frena Mt)dt Oy (...)|?dr'dy' 
“aa 04) fontane 


1 rar 
=c——___ Oy:(...)|?dx'dy', 
“nba — 6;) hi y'( )l Tay 


by (3.1), 
così che 
1 3 
i Seli -YV #(d)0:t,; 
(4 10) Ir < abi = 5)! 2,40 )p lindo» 


dove Q1, 02,3 sono i vertici di K,e | lat) denota la seminorma data dalla somma 
delle norme L? delle derivate massime. Possiamo ora applicare la stima classica 


3 
(4.11) li “ * UQ)Bil 1) < cli) 
È i=1 
(si veda [QV], Teorema 3.4.1) ottenendo 


1 2 212 222 2 112 
Wp_dai as n =; n} w 10: = 
D< r(6i =) {02 Ul) + ld, UllLa i) + ly d allz#)} 


_ ST 32,,]12 
- cf ul tati silla 


ii + MB NEU) 
= c{(6;41 = 5;) |] P a logullî:0o # =z10,2cul{t2(x} 
={7+T+J3}. 
Notiamo ora che J, può essere scritto come segue: 
Jr = (641- GY |Xulltz0): 
il prossimo passo consisterà nel provare che 
(4.13) Ta £ (641 = GEE) e Ta S e(6;4a — 6A Uli) 
A tal fine osserviamo che, se r € K, per la monotonia di À, 


1 i 1 
Messe — [ima 
26/0 


j-1 
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così che (4.13) seguirà provando che 
(4.14) 8; — 9j-1 £ (6;+1 — 8;) 


per ogni j= 0,1,...,n-1. 
D'altro canto è facile vedere che 


d((6;-1,k/m),(6;,k/n)) < 6; — 6j-1, 
così che, per la Proposizione 2.2 (ix), abbiamo: 


F((6;,k/m),5; — 6j-1) £ cF((6;-1,k/m), 6; — 6j_1) 


al i *atoa=£=o J "ma 
8j-a a 5; 

= cF((6;,k/n),6;+1 — 8;) £ F((6;,k/n), (8;41— 5;))) 
per la Proposizione 2.2 (v), e quindi (4.14) segue dalla monotonia della funzione è + 
F((6;,k/n),t). 


Dunque 
Ia < (611 — 5 {Nulla + 1248200) + RX lt } 


< n70+Df [x2ullî:co + |Xzullt20n + |xX2X1ullîa40}, 


per il Teorema 3.1 (iv) e (v). Finalmente, il termine 
h= f lu- tn(u)l'dedy 
K 


può essere stimato nello stesso modo. 


Combinando ora (4.4) e (4.5), la prova del Teorema 4.1 sarà completata provando la 
seguente stima a priori: 


— {°) 
Theorem 4.3. Supponiamo f € L°°(Q); allora esiste una ed una sola u € Hi(9) tale 
Lou=f inQ. 
Inoltre, u € C°(Q) per un opportuno o € (0,1), X?u, Xfu, X2X10 € L?(9) ed esiste 
una costante geometrica B4q > 0 tale che 
ullz + Valla + IX ullz:®) 
+ |XZu]|z2(0) + |X2X1 vl) £ Ballfilzza) 


Proof. La prova è piuttosto macchinosa, e richiede una preliminare approssimazione del 
problema con problemi ellittici in regioni regolari in modo da poter lavorare con funzioni 
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almeno tre volte derivabili con continuità. Successivamente, la prova consiste di un certo 
numero di integrazioni per parti e di argomenti ad hoc per la struttura dell’equazione 
e dell’insieme. La regolarità vicino agli spigoli segue da un risultato di Grisvard ([G], 
Teorema 1). 


5. In questo paragrafo descriveremo alcuni test numerici dei nostri risultati; come 
abbiamo sottolineato nell’Introduzione, il numero y + 2 gioca il ruolo di una dimensione, 
che può essere molto grande se l’operatore è fortemente degenere. A causa di questo, 
per stimare l'andamento delle nostre stime dobbiamo lavorare con reticoli contenenti una 
grande numero di punti, e, per questa ragione, dobbiamo scegliere accuratamente la nostra 
implementazione. L'algoritmo che abbiamo usato per effettuare la nostra integrazione 
numerica è un unified multilevel adaptive refinement method, la cui descrizione tecnica 
dettagliata può essere trovata in [M1,M7]. 

A differenza del caso ellittico, nel nostro caso il raffinamento non può essere ovviamente 
effettuato bisecando coppie di triangoli, in quanto vogliamo riprodurre la geometria associ- 
ata all’operatore differenziale considerato. Ciononostante, lo schema della base gerarchica 
può ancora essere applicato poichè, da (3.2), si può vedere che ad ogni livello di raffi- 
namento, tutti i vecchi nodi sono sono conservati, e altri nuovi sono aggiunti. Questo è 
esattamente il principio su cui si fonda l’approccio delle basi gerarchiche. Infatti la base 
gerarchica coincide con l’usuale base nodale al primo livello di raffinamento, mentre, con il 
procedere del raffinamento, uno o più nodi vengono aggiunti, e per ogni nodo viene definita 
una nuova funzione che vale 1 al nuovo nodo e 0 in tutti gli altri nodi, ma le funzioni delle 
basi esistenti non vengono cambiate. La scelta della base gerarchica porta a una rapp- 
resentazione di IIn(u) che differisce da (3.8) e quindi a una nuova matrice di stiffness: il 
guadagno algoritmico è descritto in [BDY] e [M7]. 

Come esempio, scegliamo 


\= VBmJz|® per 9>0em>1; 
si verifica facilmente che (2.4) vale con y=m — 1. Poniamo ora 
R=Bjx|" + y? e u= R!/m, 


e, per finire, 
v=(1-22)(1-y?)u. 


Un calcolo diretto mostra che v E Hi (9), e che f = Lv E L°(Q); d'altra parte notiamo 


esplicitamente che, se m > 2.5, allora v non appartiene allo spazio di Sobolev usuale H1(Q). 
Quest'ultimo fatto non é sorprendente, poiché la funzione u é equivalente al quadrato della 
distanza d del punto (x,y) dall’origine (si veda [FL]), cosi che riflette in modo piuttosto 
sottile le proprietà dell’operatore L. 

Possiamo ora valutare l’errore di discretizzazione nella norma dell’energia 


Oz(v — on)iIzza) + Ady(0 — vn)llz:(®), 
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che non solo é la sola norma associata naturalmente all'operatore, ma é anche la sola norma 
‘ragionevole’, poiché in generale la norma H! di v é infinita. Per il Teorema 2.4, la norma 


dell’energia é equivalente alla norma Hi(9). 

Nei disegni seguenti abbiamo plottato in una scala log - log questi errori come funzione di 
N = n? = numero dei nodi della nostra triangolazione (chiamata triangolazione naturale), e 
l'abbiamo confrontato con l'errore — a parità di numero di nodi — per una versione adattiva 
della triangolazione naturale, per una triangolazione uniforme e per una triangolazione 
uniforme adattiva. I grafici corrispondono am =4e m = 6. 

Abbiamo inoltre valutato il trend dell’errore confrontandolo con il valore dato dalla 
stima teorica per m = 2 e m = 3, sempre in scala log-log. In questi casi l'errore atteso é 
rispettivamente 0.25 e 0.17, che é rispettato con sufficiente precisione. 


nat en err 


uni en err 


ada en err 


adan en err 
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